On the origin and extent of mechanical variation among cells 
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Investigations of natural variation in cell mechanics within a cell population are essential to understand the 
stochastic nature of soft-network deformation. Striking commonalities have been found concerning the average 
values and distribution of rheological parameters of cells: first, attached and suspended cells exhibit power-law 
rheological behavior; second, cell stiffness is distributed log-normally. A predictive connection between these 
two near-universal findings has not been reported, to our knowledge. Here we postulate, based on our own and 
others' experimental reports and leading models of cell rheology, that the exponent that characterizes power-law 
rheology varies intrinsically among cells as an approximately Gaussian-distributed variable. Besides explaining 
naturally the log-normal distribution of cell stiffness that is widely observed, this postulate predicts multiple 
empirically observed relationships from cell deformation studies. Our framework ultimately links inherent 
noise in postulated relaxation mechanisms of cytoskeletal networks to mechanical variation among cells and 
cell populations. 
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In light of the complexity and structural heterogeneity of bi- 
ological cells, findings of universal mechanical tendencies are 
significant. Thus, it is intriguing that, on physiologically rele- 
vant time scales (at least three orders of magnitude centered on 
1 s), and in experiments at multiple length scales (Fig. [IJa)), 
animal cells regularly exhibit so-called power-law rheology 
(Fig.[TJb)). For example, dynamic cell stiffness (storage and 
loss moduli G'(qj) and G"(qj)) scales with frequency oj as 
G\<j)\G"(<jj) oc of (TH1, and creep compliance J(t) 013- 
7 ] and stress relaxation modulus G(t) 1 8 ] scale with time t as 
J it), l/G(t) oc t a , with a « 0.1-0.3 typically. This behavior is 
attributed to the varied possible arrangements in cytoskeletal 
networks, each with its own relaxation time, that integrate to 
form a material with no single characteristic time scale l9l . 

Common findings also exist regarding the distribution of 
individual cell stiffness values around a population average. 
Cell stiffness is distributed log-normally — again largely inde- 
pendent of experimental technique and length scale ITl [T0lfT2ll 
and also independent of metabolic state ifTTTl and cytoskele- 
tal perturbation (2| |3] [13) (Fig. [TJc)). Wide population dis- 
tributions of stiffness measurements have previously been at- 
tributed to probe contact variation QEL but arise even with 
non-contact single-cell deformation 0. In contrast to the ex- 
istence of multiple models of power-law rheology (PLR, re- 
viewed in [14]), the observed log-normal distribution has re- 
mained largely uninvestigated and especially lacks any pre- 
dictive explanation. (A previous report replicates the distribu- 
tion shape via a random selection of viscoelastic elements, but 
does not provide a way to predict distribution parameters l5l.) 
Also unexplained are observed differences in the width of dis- 
tributions (quantified by geometric standard deviation SD geo ) 
between storage and loss modulus measurements lfT2l . As we 
develop a better understanding of the cell as a mechanical ma- 
terial, natural questions are (1) whether PLR implies a log- 
normal distribution (and if so, under what necessary postu- 
lates); (2) what is the minimum measurable mechanical vari- 



ation in a cell population; and (3) whether fluctuations in the 
energy landscapes of cytoskeletal networks can be linked to 
mechanical variation among cells. 

With an interest in answering these questions, we propose 
that the power-law exponent a varies intrinsically among cells 
approximately as a Gaussian-distributed variable, based on 
the following evidence: First, experimental measurements of 
the exponent a have actually been found to exhibit a Gaus- 
sian distribution (4l |5] HQj [T21. Moreover, literature reports 
on the variation of a (quantified as standard deviation cr a ), as 
shown in Fig. [2j suggest an intrinsic variation that lies in the 
approximate range 0.02-0.10 and is sustained even when hun- 
dreds of individual cells are sampled. Second, multiple PLR- 
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FIG. 1 : (a) Universal mechanical behavior of living cells at multi- 
ple length scales: (b) The storage and loss moduli G'{co) and G"(oj), 
along with the creep compliance J(t) and reciprocal stress relaxation 
modulus 1 /G(f), scale as a power law with frequency a> or time t, 
respectively; (c) Individual measurements of these rheological pa- 
rameters are distributed log-normally. 
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FIG. 2: Existing estimates of standard deviation cr a of power-law 
exponent a as a function of data set size suggest that the inherent 
uncertainty in a lies between 0.02-0.10. Data referenced by study 
author and year t3l[5ll6lfTTWT3l[T7l. 



predicting models of mechanical networks relate a to an en- 
ergy. (For example, in the soft glassy rheology (SGR) model 
developed by Sollich et al. l2T1l and applied to cells by Fabry 
and Fredberg et al. [Q] 0, x = a + 1 is the ratio of the av- 
erage agitation energy to the energy associated with a glass 
transition; in the glassy wormlike chain model, S ~ 3/a is 
the average energy barrier retarding material relaxation l22l .) 
Notably, these models assume that the power-law exponent 
represents the sum of many independent energies (21] [22) ; 
from the central limit theorem, therefore, we would expect 
the exponent a to be approximately Gaussian regardless of 
the distributions of its constitutive energy components. Thus, 
our postulate of intrinsic Gaussian exponent appears plausi- 
ble in both experimental and theoretical contexts. We fur- 
ther assume that endogenous variation in a dominates over 
other sources of variation such as direct engagement of the 
cytoskeleton via adhesion complexes 0. 

We now relate postulated fluctuations in the power-law ex- 
ponent a to quantitative predictions of variation in mechanical 
parameters such as the complex modulus G*(oj) = G'(qj) + 
iG"(o)), creep compliance J(t), and stress relaxation modu- 
lus G(t). Consider first the complex modulus and its stor- 
age and loss components, expressed as f(a)(oj/Yo) a , where 
f(a) is a characteristic prefactor for that component (namely, 
T(l+a)T(l-a) cos[7r(a/2] for G' mdT(l+a)T(l-a) sin[7ia/2] 
for G") as obtained from the structural damping and SGR 
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FIG. 3: In non-contact measurements of individual mesenchymal 
stem cells by optical stretching, stiffness G is negatively correlated 
with power-law exponent a for sinusoidal measurements at 0.5 Hz. 



models |20l|2Tj, and where Yq is a scaling frequency. In the 
SGR model, for example, Yq represents the maximum possi- 
ble yielding frequency of a mesoscopic region I2l1 . We as- 
sume that a is Gaussian with average a and standard deviation 
a a for a given study (see Fig. [2] for reported studies). We use 
the change-of-variables equation P[G(a>)] = \da/dG(a>)\P(a) 
to obtain 



P[G(oj)] = 



1 



1 



GM^ L (in 



Yo _ dlnfja) 
co da 



exp 



[lnG(oj)/f(a)-lnG(oj)/f(a)] 



(1) 



a variant of the log-normal distribution in which the place of 
the standard deviation is taken by the absolute value term. 
When f(a) is only weakly dependent on a, as is the case 
for the storage modulus (Table [I]), we have dlnf(a)/da « 
and the standard deviation is approximately cr a ln(Yo/cj). 
This standard deviation corresponds to a geometric standard 
deviation SD geo « (Yo/oj) (Ta for this mechanical parameter. 
When this approximation is not valid (specifically, for the 
loss modulus), and when din f (a) /da < ln(Yo/a>) (which ap- 
plies for typical experimental condition oj Yq), the base 
term (Yq/oj) ' is multiplied by a sub-unity correction factor 
exp[-cr a din f (a) /da] « exp(cr a /a). We list the resulting 
complete SD geo terms and first-order, Taylor- series-expanded 
approximations for a > in Table [I] 

We support the model's framework with experimental data 
obtained by optical stretching of cells in the fully suspended 
state; this technique involves no physical probe-cell contact, 
making it a favorable method for evaluating intrinsic cell-to- 
cell mechanical variation without conflating variation from 
probe contact. Adult human mesenchymal stem cells were 
cultured and deformed, and cell deformation analyzed, as de- 
scribed previously (7), except that the cells were interrogated 
in the frequency domain by irradiating them with a 0.5 Hz si- 
nusoid with a mean power of 1 W per fiber and a peak-to- 
peak value of 1 W per fiber. Cell stiffness G was calculated as 
ctq/sq, where cr is the amplitude of the photonic stress and 
sq is the amplitude of a sinusoid fitted to cell deformation; the 
power-law exponent a was calculated as 2(p/n where cp is the 
phase lag of cell deformation. 

As shown in Fig. [3] our oscillatory measurements produced 
a trend of decreasing stiffness G with increasing power-law 
exponent a. This relationship is predicted by SGR theory 
(because a larger agitation energy implies relatively more de- 
formation for a given load [21]). Such a prediction has not 
been demonstrated previously for a single experimental con- 
dition due to the conflation of probe attachment variation 
that is possible with other approaches; in fact, this probe en- 
gagement has been reported to introduce the opposite corre- 
lation O [13] |20|. Balland et al.'s model of log-normal dis- 
tribution origin, which employs randomly selected groups of 
viscoelastic components, replicates this opposite correlation 
(i.e., increasing G with increasing a) |5|. In contract, our fo- 
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TABLE I: Predicted geometric standard deviation SD geo of log-normal distributions for different mechanical parameters in the context of cell 
power-law rheology, where complex modulus increases with frequency as of and creep compliance and reciprocal relaxation modulus increase 
with time as t°, where the power-law exponent a > 0. 
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cus in this work is inherent cell-to-cell mechanical variation 
alone; in this context, the trend we show in Fig.[3]and describe 
mathematically above is internally consistent and is also com- 
patible with the SGR framework. 

Having constructed a connection between the power-law- 
uncertainty cr a and the geometric standard deviation SD geo 
of mechanical measurements for complex modulus G*(cS) = 
G'(qj) + iG"(aj), we now compare our resulting predictions to 
experimental findings from other groups, represented in Fig. [2] 
and [4] Four predictions follow: First, a log-log plot of SD geo 
vs. measurement period 1/cj will have a slope cr a for any 
group of measurements. Second, extrapolation to SD geo = 1 
on this plot will occur at time I/Yq. Third, the correction fac- 
tor derived above will be approximately 1 for G' and approxi- 
mately exp(-cr a /a) for G" (at frequencies oj <?c Y , where this 
restriction is explained in the derivation above). For example, 
we predict that the SD geo for G" will be approximately 20% 
less than for G' for a « 0.20 and cr a « 0.05. (The lower vari- 
ation of G" relative to G' can be understood by considering 
the hysteresivity G n '/G' = tan(7ra/2) in the structural damp- 
ing model. Larger values of a reduce stiffness, but increase 
G" relative to G f . Therefore, fluctuations in a are naturally 
suppressed in G".) Fourth, SD geo will generally increase with 
increasing period 1 /<d; conversely, the distribution of mechan- 
ical measurements will become more narrow with increasing 
measurement frequency. 

All four of these predictions are confirmed by results and 
trends previously reported in the literature but unexplained up 
to this point (Fig. [2] and [4]). We have extracted cr a and SD geo 
values from multiple cell rheology reports, including a set of 
oscillatory indentation measurements of G' and G" by Hirat- 
suka et al. fT2ll . Figure [4] shows a plot of SD geo values of G' 
and G" extracted from Hiratsuka et al.'s presentation of dis- 
tributions at different frequencies, along with other groups' 
values. By fitting our model, we obtain the slope and unity 
intercept estimates cr a « 0.05 and 1/Fq ~ 1 ills. The first fitted 



parameter is in good agreement with our estimate for intrin- 
sic exponent variation described above, with cr a slightly less 
than the standard deviation (0.064, see Fig. [2]) of that group's 
reported a values. This finding supports our first prediction 
above. Second, the estimate of Y$ « 1 MHz is in fairly good 
agreement with Fabry et al.'s estimate l25ll . Third, Hiratsuka 
et al. found SD geo for G" to be 14-25% less than that for 
G r (Fig.|4]), also in agreement with our prediction of how SD geo 
should vary between the real and imaginary components of the 
complex modulus G*. Fourth, the reported or extracted SD geo 
generally does increase with period for SD geo values reported 
in the literature. 
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FIG. 4: Estimates of intrinsic variation lead to predictive models. 
The observed variation in cell stiffness moduli, as quantified by geo- 
metric standard deviation SD geo of a log-normal fit, as a function of 
experimental time scale, for studies of >100 cells. Solid lines show 
our predictions relating SD geo to frequency for storage (upper) and 
loss (lower) moduli for measurements of the same cell population 
at different frequencies fT2l . according to equations described in the 
text. The slopes correspond to cr a = 0.05 (best fit), in agreement 
with the estimate made in Figj2 Data referenced by study author 
and year El [6] [3 El [Bl [T6lll8ll201; bolded references l71[T2l[T9l 
show studies not requiring transmembrane probe-cytoskeleton link- 
ages. (The stochastic nature of cytoskeletal engagement and anchor- 
ing tends to introduce additional experimental variation \13\ [20) 
not addressed by our model.) 
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The distribution widths for creep compliance J(t) and re- 
laxation modulus G(t) are readily calculated in the same way 
as the complex modulus. Different pref actors are possible for 
each of J(t) and Git) depending on whether t «: t w or t » t w , 
where t w is the time between a fluidizing large-strain event and 
the start of the experiment, as modeled by Fielding et al. l23l 
and experimentally explored in live-cell studies conducted by 
Bursac et al. and Trepat et al. [021 [24). Despite these addi- 
tional scenarios, the determination of SD geo parallels our ear- 
lier treatment of the complex modulus. The resulting SD geo 
terms of G(t) and J(t) and first-order, Taylor-series-expanded 
approximations for a > are listed in Table [I] as predictions to 
be verified [26]. 

Revisiting the natural questions arising from universal ob- 
servations of cell populations, we find that an assumed dis- 
tribution of power-law exponent a — together with prefactors 
calculated from phenomenological and theoretical models of 
cell rheology — lets us calculate the distribution of dynamic 
stiffness, creep compliance, and stress relaxation modulus. A 
Gaussian distribution of a leads immediately to a log-normal 
distribution of stiffness (and of compliance, due to the recip- 
rocal nature of this distribution). We also find the minimum 
measurable geometric standard deviation SD geo for experi- 
ments conducted in the frequency regime to be oc(Yo/a))°~ a . 
The prefactor here depends upon the experimental regime, 
mechanical parameter of interest, and cell state, as listed in Ta- 
ble[l| We have presented a new way to estimate cr a and Yq from 
the log-log slope and unity intercept of SD geo as a function of 
experiment frequency. (Note that fuller interpretation of these 
estimates is enabled when SD geo is reported over a range of 
frequencies in a given study, as in Ref. 12.) Additionally, we 
have predicted and confirmed an unexpected relationship be- 
tween the distribution widths of G' and G" measurements for 
a given cell population. Finally, given that the parameter a 
has been linked to average activation energy or barrier height 
in the energy landscapes of so-called soft glassy materials, we 
conclude that fluctuations in a (among cells and/or over time 
in individual cells) could plausibly produce the cell-to-cell 
mechanical variation that is universally observed. Our analyt- 
ical model differs from previous investigations [ 5 ] that repro- 
duce distribution shapes via phenomenological viscoelastic 
models, now providing testable predictions that relate phys- 
ical mechanisms (e.g., attempt frequency) to measured distri- 
butions. 

Our framework is applicable to soft glassy regions — 
animate and inanimate — that exhibit PLR, as long as the stan- 
dard deviation of the power-law exponent quantifies the dom- 
inant source of mechanical variation. These requirements are 
concluded herein to be met by living cells that are measured 
individually and reported as a population distribution that is 
consequently log-normal. Our model focuses on endogenous 
variation and its manifestation as measured by low-contact 
or no-contact techniques (such as oscillatory probe indenta- 
tion [12] and optical stretching [7], respectively), in which 
probe-cytoskeleton linkages are not required, thus reducing or 
eliminating measurement variation caused by ligand-receptor 



engagement of the cytoskeleton. Additional experiments and 
models are needed to determine conclusively whether disease 
state or chemomechanical modulation of the cell environment, 
both of which are known to alter cell stiffness, can also modu- 
late the mechanical variation among cells [ 27 ]. Most valuable 
are high-throughput cell measurement techniques that allow 
robust estimates of distribution width as a function of time 
or oscillation frequency, as well as estimated errors in these 
widths. These quantities will inform models such as ours that 
can then connect microscopic structural and energetic barriers 
within cells to emergent changes in mechanical distributions 
among cells. 
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